Detecting charge noise with a Josephson junction: 
A problem of thermal escape in presence of non-Gaussian fluctuations 
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Motivated by several experimental activities to detect charge noise produced by a mesoscopic 
conductor with a Josephson junction as on-chip detector, the switching rate out of its zero-voltage 
state is studied. This process is related to the problem of thermal escape in presence of non-Gaussian 
fluctuations. In the relevant case of weak higher than second order cumulants, an effective Fokker- 
Planck equation is derived, which is then used to obtain an explicit expression for the escape rate. 
Specific results for the rate asymmetry due to the third moment of current noise allow to analyse 
experimental data and to optimize detection circuits. 

PACS numbers: 72.70.+m,05.40.-a,74.50.+r,73.23.-b 
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A complete understanding of electronic transport 
through mesoscopic conductors necessitates the knowl- 
edge of all noise properties of the corresponding current. 
This is the goal of full counting statistics, which has at- 
tracted considerable activities in the last years [j], 
Theoretically, generating functionals of current noise in 
the low frequency limit have been calculated for a vari- 
ety of systems comprising tunnel contacts, diffusive wires, 
and ballistic cavities |2|. Experimentally, in a pioneering 
measurement the third moment of current noise produced 
in a tunnel junction was detected by analog amplifiers 
and filtering techniques in Q and later also in Since 
then strong efforts have been made towards on-chip de- 
tection schemes, first because they are faster and second 
because they give access to finite frequency noise proper- 
ties. Lately, the distribution of charges flowing through 
a quantum dot has been extracted on-chip in the low fre- 
quency regime The goal now is to push devices into 
higher frequency ranges (GHz), where quantum effects, 
electron-electron interactions, and plasmon dynamics are 
relevant. 

Based on an idea proposed in [6| and studied later in 
various scenarios 0, [1[ , currently, several experiments are 
aiming to set up circuits with Josephson junctions (JJs) 
as detection elements [H, [l(| EH • JJs can be fabricated 
and manipulated in a very controlled way, offer a large 
bandwidth depending on their plasma frequency, and 
contain an intrinsic amplification mechanism. Namely, 
the switching out of the zero voltage state is exponen- 
tially sensitive to variations of the barrier potential and 
to the noise strength. For electrical noise large positive 
and negative fluctuations from the mean occur with dif- 
ferent probabilities. The idea is thus to probe this asym- 
metry by measuring switching rates for mean mesoscopic 
currents flowing forward and backward, respectively. In- 
deed, recent experimental results 111] indicate a sufficient 
sensitivity of such a circuit to retrieve the third cumulant. 

A typical set-up [l(| consists of a J J, on which two 
currents lb and I m are injected. Current lb is a stan- 
dard bias current coming from a source in parallel to 
the J J, while I m runs through a noise generating meso- 
scopic conductor in series with the JJ in such a way 



that no dc component of I m passes through the detec- 
tor. Due to substantial heating from the additional elec- 
trical noise the JJ operates in the regime of classical es- 
cape, where quantum effects are negligible. Then, ac- 
cording to the resistively and capacitively shunted junc- 
tion model, the phase tp of the JJ moves in a tilted wash- 
board potential —Ej cos (ip) — (h/2e)(Ib)<p with Joseph- 
son energy Ej and is subject to Johnson-Nyquist noise 
51b = lb — (h) and stationary non-Gaussian current fluc- 
tuations SI m = I m —{I m ) . If the phase is initially trapped 
in one of the wells (zero- voltage state) , it may for suffi- 
ciently large (h) < I c = (2e/fr)Ej escape so that the J J 
switches to a finite voltage state. Further, since the third 
cumulant vanishes in equilibrium due to time-reversal 
symmetry, experimentally, the mesoscopic conductor is 
at low temperatures driven far from equilibrium into the 
shot noise regime, where no fluctuation-dissipation theo- 
rem applies. Hence, the switching of the JJ can be visu- 
alized as the diffusive dynamics of a fictitious particle in 
a metastable well with non-Gaussian continuous fluctua- 
tions acting as external random driving force. 

In the past JJs have been used to thoroughly con- 
firm various theories for thermal escape [l2l [l3j including 
cases of external driving either by time-periodic forces 
or Gaussian noise sources related to phenomena such 
as resonant activation [3, EBl and stochastic resonance 
[l6| . When putting the experimental situation described 
above into this context, one arrives at a new type of rate 
problem: Thermal escape driven by non-Gaussian con- 
tinuous noise. Theoretically, the challenge here is that 
externally driven escape necessitates a full dynamical de- 
scription because the stationary well state from which 
particles are ejected over the barrier is only known a 
posteriori. In contrast, in the undriven situation ther- 
modynamic approaches, as e.g. the bounce method, can 
rely on a thermal equilibrium state inside the well region. 
In the classical domain and for Gaussian noise processes a 
dynamical formulation is based on a Fokker-Planck equa- 
tion (FPE) for the phas e sp ace distribution as shown in a 
variety of applications [LJ, |l7|. Thus, the first goal is to 
derive a generalized FPE for weak non-Gaussian noise, 
which then serves as a basis for the rate calculation. A 
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solution to this problem is of general interest for rate pro- 
cesses in complex media and crucial for present electrical 
noise measurements because non-Gaussian components 
are typically very small compared to a prevailing Gaus- 
sian background. In this Letter we develop the frame- 
work for such a rate theory and give explicit expressions 
in case of a JJ as detector for electrical noise. 

For this purpose let us consider the Langevin equation 
for the diffusive motion of a particle of mass m 

mip{t) + V\ip) - v(t) + m 7 p{t) = , (1) 

where ip denotes a generalized coordinate, ' = d/dt, and 
' = d/dip. The barrier potential V(<p) is assumed to 
be sufficiently smooth with a we ll located at ip = 
with frequency ujq = y/V"(0)/m and a barrier top at 

if = ifb with frequency u>b — \Z\V'( i Pb)\/i r n and height 
Vb = V(ifb) — V(0). The white thermal Gaussian noise 
obeys (£(t))p = and is related to the damping 7 
via the dissipation fluctuation theorem: (£(t)£(t'))/3 = 
(27717//?) S(t — t') with inverse temperature /? = l/k^T. 
The statistical properties of the stationary non-Gaussian 
noise are determined by the generating functional 
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where T is the time ordering operator and the cumulants 
are gained from the functional derivatives of S[w] . In par- 
ticular, we assume Si(t) = (r)(t)) — i dS[w]/dw(t)\ w=0 = 
0, while the auto-correlation function reads 

s 2 (t) = ( v (t)v(o)) = d 2 s[w]/dw(t)dw(o)\ w=Q 

and the third cumulant follows accordingly as 

S 3 (t,t') = {r 1 (t + t')r 1 (t , )r ] (0)). 

Now, particles initially confined in the well region may 
escape such that for sufficiently high barriers (/3V& 3> 1) 
a stationary flux appears related to an escape rate T — 
(v5(<p — (/5b))flux/^Vwcii- The denominator is the popula- 
tion in the well and the expectation value of v = ip is 
taken with respect to a quasi-stationary noncquilibrium 
state, the flux state, which may be cast into the form 



L (v,ip) = P cq ui(v, ip) /fl ux (w, <p) 



(3) 



Here P e qui denotes the equilibrium distribution and /a ux 
describes deviations from it such that it tends to to- 
wards the continuum and to 1 towards the well. The flux 
state, as a phase-space density, is obtained as a station- 
ary solution to the Fokkcr-Planck equation (FPE) corre- 
sponding to ([T]). In case of vanishing non-Gaussian noise, 
?y(i) = 0, the latter one reads dtP(v,<p,i) = LqP(v, p,t) 
where 

L a = -vd v + d v [V'(<p)/m + jv] + 7 /(™/?R 2 . (4) 

Then, from Lo^flux = one finds for moderate 
to strong friction the known expression [13] T = 



exp(— (3Vb) with the scaled Grote-Hynes fre- 
quency n = -7/(2w 6 ) + y/(j/2u b ) 2 + 1. 

For finite non-Gaussian noise the translation of |T]) into 
an equivalent FPE for the averaged phase s pac e distri- 
bution leads to only formal expressions (l7l. Il8j|. Here, 
in accordance with the experimental situation we pro- 
ceed by assuming that (i) non-Gaussian fluctuations are 
weak and (ii) sufficiently fast compared to the bare dy- 
namics of the system (detailed conditions will be given 
below). Then, one considers V e s = V — iprj(i) as an effec- 
tive time-dependent potential so that for each realization 
of the non-Gaussian random force a FPE for P v (v,ip,t) 
exists where Lq is replaced by Lq — T](t)/md v . To gain a 
FPE for the averaged distribution (P v (v, ip, t)), we switch 
to the interaction picture P v (t) = exp[Lo£] Q n (t) with 
d v (t) = exp(— Lot)d v exp(Lot). By averaging the equa- 
tion for Qr)(t) over the 77-noise, the exact solution is ex- 
pressed in terms of the generating functional with the 
counting field w(s) substituted by (i/m)d v (s), i.e., 



(Q„(t)) = c 

= Texp 



S[(i/m)d v (s)] q(q) 

C k (s) 



J ° k>2 



0(0). (5) 



Here, we have used the cumulant expansion of 
S[(i/m)d v (s)] with Cfc(s) = exp[—L s]Ck(s)exp[L s], 
where the two lowest order cumulant operators read 



C 2 (s) 



C 3 (s) 



duS2(u)d v d v (—u) 

S pS~U 

du I du' Sz(u,u')d v d v {—u) d v (—u — v!) 
<< Jo 



The above result reveals that the generator for the aver- 
aged dynamics is directly given by the generating func- 
tional of the non-Gaussian noise. Its cumulant expansion 
now allows for a systematic approximation. 

To do so, we take only the two leading order terms 
into account in which in turn leads to d t (P v ) — 

L + C 2 {t)/m 2 + C 3 (t)/m 3 (P v ). While this FPE ap- 
plies to weak non- Gaussian noise of arbitrary type, fol- 
lowing assumption (ii), the noise correlation functions 
vanish on a sufficiently short time scale so that the non- 
Markovian dynamics reduces effectively to a Markovian 
one. This way, one arrives at a generalized FPE of the 
form d t (P) =L eS (P) with 



i cff = L Q + (c 2 /m 2 )d 2 v - (c 3 /m 3 R 3 



(0) 



where 



c 2 = I duS 2 (u) , C3 = / du du' S 3 (u,u') . 
Iq Jo Jo 

Note that this expression generalizes previous results ob- 
tained in the overdamped limit [l8j to weaker friction. If 
the external noise were purely Gaussian, i.e. C3 = 0, the 
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operator L e g were exact. Hence, what we really have to 
assume in deriving L e g are small higher than second or- 
der cumulants, while c 2 may be arbitrarily large. Thus, 
one introduces an effective temperature 

Trff = T + c 2 /(k B mj) (7) 

and incorporates the Gaussian components of the non- 
Gaussian noise into a renormalizcd diffusion term accord- 
ing to Lq{0) + (c 2 /m 2 )d 2 — > L (/3 c ff). We note in pass- 
ing that this expression coincides with the one derived 
for resonant activation [14|. Experimentally, the heating 
due to c 2 is substantial so that J7]) is required to capture 
the actual temperature of the J J. 

With the generalized FPE at hand we now attack 
the rate calculation. For this purpose, it is convenient 
to work with dimensionless quantities r = oj$t, x = 
ip^P eS mu} 2 , p = v^j3 eS m, U = f3 cS V, and p = 7/wo, 
C3 = C3(/3 e ff/"i) 3 / 2 ; then Xb,Ub ~> 1 and C3 serves as a 
small parameter. We start with the equilibrium state (no 
flux) which determines the dominating exponential acti- 
vation factor in (J3|) and write P equ i oc Pp B{! exp(— C3G) 
with the Boltzmann distribution P/3 eff . Upon inserting 
this ansatz into L c ffPcqui = we find perturbatively a 
distribution of the form 

3 

G(x,p) = <j> (x) + J2 ( l , n(x)P n /n> (8) 

71=1 

where coordinate dependent functions can be expressed 
in terms of only 4> 2 as 

Mx) = -p<h(x)/U'(x) = -fa(x) 

pX pX 

4> {x) = 3x-3p dyMv)+ dyU'(y)My)- (9) 
Jo Jo 

Note that the property 4>\ = —tfiz ensures that (p) oqu i = 
in leading order. While the above expressions hold in gen- 
eral, the solution of 4>2 depends on the specific form of the 
potential U. To better understand the deviations from 
the thermal equilibrium wc thus first study a harmonic 
potential U(x) = x 2 /2, which approximates a metastable 
potential around its well minimum. The result is 4> 2 (x) — 
—2x/(l + 2p 2 ) and the corresponding distribution P equ i 
has a simple interpretation: it coincides with the thermal 
distribution driven by the external force 77(f) and aver- 
aged in the long time limit over the non-Gaussian noise. 
Specifically, one has P cq ui(p, x) = lim T _ f00 (P / 3 etf (p T , x T )), 
where p T = p — ,/(r) and x T — x — /(r) with /(r) — 

\J /? c ff / J °° rfr'x(r — t')t](t') and x( T ) the response 
function for a damped harmonic oscillator. As a conse- 
quence, P C qui {p, x) shows an asymmetry even for a sym- 
metric potential depending linearly on the third cumu- 
lant S3. In contrast, in a set-up where for the r\ noise 
a dissipation fluctuation theorem and correspondingly a 
detailed balance relation applies, a linear dependence on 
odd cumulants only occurs on a transient time scale 

A generic metastable potential is of the form U (x) = 
(x 2 /2)(l — 2x/3xb) describing particularly a well-barrier 



segment of a tilted washboard potential of a J J. Then, 

, . . 2xhz(l — z) a , . , . 

<t>i{x) = 1 — 2 Fi(a,l + a,2 + a,z) 10 

1 + a 

where z = x/xb, a = 2p 2 and P equ i follows together with 
([9]). There is a little subtlety here in that the solution 
(HHJ) applies for 2p 2 < 1 only outside a narrow range 
around Xb- The global solution is then constructed by 
properly matching the local solution around Xb onto the 
latter one. Since the vicinity around the top affects only 
the prefactor and not the experimentally dominant ex- 
ponential factor of the rate, we will give further details 
elsewhere. This exponential activation factor can now 
be inferred from P oq ui(p = 0, x = Xb) and is obtained as 
r oc exp[— Ub(l — g)] with the third cumulant correction 
g = — cz4>o{xb) /Ub explicitly evaluated as 

The remaining prefactor to the rate is determined by 
/flux in ([3]) and the well population N wc \\. The latter 
one is easily calculated from P equ i(j3, x) restricted to the 
domain around the well bottom, while the former one is 
derived from the local dynamics around the barrier top. 
Eventually, the rate is found as 

where k = 20(9 + 4p 2 ) + p[8p 2 + ^/p^{7 + 2p 2 )} + 4 and 
4>' 2 (xb) results for 2p 2 < 1 from the matching procedure 
leading, e.g., for 2p 2 <C 1 to <j)' 2 {xb) = 1/p 2 - The above 
expression is the main finding of this Letter, namely, the 
thermal escape rate out of a metastable well in presence 
of weak (03 <C l/£4) and fast (correlation time r c -c 
1/wo) non-Gaussian continuous noise. This result indeed 
verifies that the detector transforms the noise asymmetry 
into a rate asymmetry depending on the sign of the third 
cumulant C3. For C3 > 0[< 0] the noise distribution favors 
large positive (negative) fluctuations so that the barrier 
is effectively reduced (enhanced) and the distribution in 
the metastable potential develops a tail towards large 
positive (negative) momenta and coordinates causing a 
rate increase (decrease). 

In the sequel, the general result (L2"|) is applied to a cir- 
cuit, where a JJ acts as detector for the current noise pro- 
duced by a normal tunnel junction subject to a voltage V 
in the shot noise regime eV 3> k^T. Further, eV/% » ujq 
guarantees that the noise is much faster than the plasma 
frequency uj = (y/2EjE c /K)(l - s 2 ) 1 / 4 (s b = h/h) of 
the JJ with charging energy Eq = 2e 2 /C (capacitance C) 
and that its back-action onto the tunnel junction is neg- 
ligible. Then, the current statistics is purely Poissonian: 
during a (scaled) time interval t p an average number of 
N charges passes the conductor so that s m = (I m /I c ) = 
Neuj /{T P I C ) and S 3 (t,t') = 5{t + r')<5(r)(ewo// c ) 2 s m . 
The procedure to extract the third cumulant is to mea- 
sure switching rates with bias currents Sb (or equivalently 
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FIG. 1: Rate asymmetry Rr according to (fl3|) lot luqTiP = 0.5, 
/3Ej — 200, Q — 2.5, and various s m vs. the bias current st,: 
s ra = 3 (solid), s m ~ 2.3 (long-dashed), s m = 1.7 (short- 
dashed), s m = 1.1 (dotted); Inset shows Rr for the same 
parameters and Sb = 0.75, s m — 1.1 as a function of Q. 

mesoscopic currents s m ) flowing forward and backwards, 
respectively fllj . This way, one finds for the rate asym- 
metry i?r = r(|s;,|)/r(— \sb\) as the dominant contribu- 
tion Rr as exp[2C/fc<?(|sh|)]. In the case considered here 
the result for this asymmetry expressed in junction pa- 
rameters reads 



Rr = exp 



AV2(3^EcE 2 jQ 2 a m (l- |s 6 |) 2 



3(5 + Q 2 ) 



Sb 



(13) 



l/p and the effective tem- 
l /(4k By /l - si). Note that 
in the exponent the ratio between third cumulant and 



with the quality factor Q = 
perature T e g = T + QTlluq s, 



Gaussian noise basically appears via Ej(3 2 s s m so that 
Rr — > 1 for increasing s m . The above finding not only 
allows to understand experimental data, but may be used 
to optimize the detection circuit as well. Namely, as a 
function of Q, the asymmetry Rr exhibits a maximum 
for intermediate Q-values, but decreases towards higher 
(underdamped) and lower (overdamped) ones (Fig. [T]). 
This also reveals that a rate calculation in either of these 
limiting cases is not sufficient. For typical experimental 
parameters s m = 2lu Ti(3 = 1, (3Ej = 200, and s& = 0.75 
one gains a maximal Rr ~ 1.45 at Q w 2.5. Asymmetry 
ratios for various s m , see Fig. [TJ lie around this value in 
agreement with recent experimental findings [ll|. Note 
that for these parameters T e g/T as 1.5. Precise numerical 
simulations confirm the result (fl3|) and will be discussed 
elsewhere [lj|. The general results (|11[) . (IT21 determine 
the rate asymmetry due to the third cumulant also for 
other mesoscopic conductors such as diffusive wires or 
chaotic cavities. Since they apply for fast fluctuations, 
they provide a tool to analyse current noise measure- 
ments in the interesting high frequency regime (2p| | . 

To summarize, we have developed a formalism to de- 
scribe the switching process of a JJ in presence of weak 
and fast non- Gaussian fluctuations. From the corre- 
sponding rate expression the rate asymmetry due to the 
third moment of current noise has been obtained. Our 
results for the generalized FPE ((6]) and the escape rate 
(fT2|) are applicable to other decay processes in physics 
and chemistry, where complex noise sources are present. 
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